OR-iS3l st>t 


FEASIBILITY STUDY OF 
TRANSIT PHOTON CORRELATION ANEMOMETER 
FOR AMES RESEARCH CENTER UNITARY WIND TUNNEL PLAN 


FINAL REPORT 


by 

W. T. Mayo, Jr. 
A. E. Smart 


7 February 1979 


SDL No. 79-6397 


Prepared for: 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Ames Research Center 
Moffett Eield, California 94035 

Attention: Dr. Dennis Johnson 

Under Contract No . : 

NAS2-10072 


&BI6HS& SOOTHS 
©SUE IlLuiTSAHMS 


PECTRON 

Development 

ABORATORIES 




APR1970 

SECE|y F0 


JLMC. 


NASA 


sn 


4co Cri*v 


(NAS A-CR-152238) FEASIBILITY STUDY OF 
TRANSIT PHOTON CORRELATION ANEHONETER FOR 
AMES RESEARCH CENTER UNITARY HIND TUNNEL 
PLAN Final Report (Spectron Development 
Labs., Inc.) 86 p BC A05/HF A01 CSCL 14B G3/09 

3303 Harbor Boulevard, Suite G-3 
Costa Mesa, California 92626 (714) 549-8477 


facilit y 


cnl 


.55 


N79- 


Unclas 

16699 




TABLE OF CONTENTS 


Page 


1.0 INTRODUCTION 1 

2.0 HISTORY AND DESCRIPTION OF SDL LTA SYSTEM 3 

2.1 Historical Introduction 3 

2.2 Description of System Used for 

NASA Ames Experiment 6 

3. 0 MEASUREMENTS AT NASA AMES 14 

3.1 Pre-Tunnel Procedures . 14 

3.2 Tunnel Measurements: Group 1, 

October 16, 1978 16 

3.3 Tunnel Measurements: Group 2, 

October 17, 1978 17 

3.4 Data Presentation: Group 2 "18 

4.0 DATA ANALYSIS AND INTERPRETATION 31 

4.1 Definition of Data Rate 32 

4.2 Data Rate and Background Level 

Equations 35 

4.3 Data Rate, Scaling with Laser Power 

and Optical Efficiency 36. 

4.4 System Particle Size Sensitivity 

Calculations ' 39 

4.5 Geometric Scaling Laws 50 


-i- 



Page 


5.0 SUMMARY AND DESIGN RE COMMENDATION S 53 

6 . 0 REFERENCES 55 


APPENDIX — PLOTS OF THE CORKELOGRAMS FROM 57 

.RUNS OCTOBER 17, 1978. 



LIST OF FIGURES 


Figure 


Page 

1 

OVERVIEW SDL LASER TRANSIT ANEMOMETER 
(LTA) OPTICAL HEAD 

7 

2 

SCHEMATIC DIAGRAM OF OPTICAL COMPONENTS 
IN SDL OPTICAL HEAD 

8 

3 

SDL LTA ELECTRONIC COMPONENTS 

9 

4 

BLOCK DIAGRAM OF SDL LTA SYSTEM FOR 
AMES TRANSONIC MEASUREMENTS 

10 

5 

CURVE FIT OF COKRELOGRAM AT M = 0.9, 
THRESHOLD = 400 mv 

19 

6 

CURVE FIT OF CORRELOGRAM AT M = 0.9, 
THRESHOLD = 280 mv 

20 

7 

CURVE FIT OF CORRELOGRAM AT M = 0.9, 
THRESHOLD = 200 mv 

21 

8 

CURVE FIT OF CORRELOGRAM AT M-= 0.9, 
THRESHOLD = 140 mv 

22 

9 

CURVE FIT OF CORRELOGRAM AT M = 0.9, 
THRESHOLD = 100 mv 

23 


-ix— 



LIST OF FIGURES 


Figure Page 


10 CURVE FIT OF CORRELOGRAM AT M = 0.9, 

THRESHOLD = 70 mv 24 

11 CURVE FIT OF CORRELOGRAM AT M = 0.9, 

THRESHOLD = 50 mv 25 

12 APPARENT TURBULENCE INTENSITY ON CONDITION 

AT M = 0.9, Re. No. - 3xl0 6 /ft 29 

13 OSCILLOSCOPE TIME EXPOSURES OF THE 
ELECTRONIC ADDITION OF CHANNELS A AND B 
PHOTOMULTIELIER TUBE OUTPUT MONITORS, 

OCTOBER 16, 1978 33 

14 DATA PAIR RATE VERSUS THRESHOLD FOR 

MACH 0.9 RUNS, OCTOBER 17, 1978 37 

15 LOG-LOG PLOTS OF DATA RATES FROM 

FIGURE 14 38‘ 

16 BACKS CATTER MIE COEFFICIENTS FOR 

_ SUBMICRON SIZED DROPLETS, n = 1.4+jO 40 

17 SIMULATED PHOTOMULTIPLIER SIGNALS, 

0,5xl0 -6 m PARTICLES, OPTICAL 

EFFICIENCY 0.3 44 

18 LOW-PASS FILTERED SIGNALS FROM 

FIGURE 17 45 

19 SIMULATED SIGNALS, 0.5xl0 -6 m PARTICLES, 

OPTICAL EFFICIENCY 0.3x0.15 (PULSES 

LOCATED AT SAME TIMES AS IN FIGURE 17) 46 

20 LOW-PASS FILTERED SIGNALS FROM FIGURE 19 47 

21 SIMULATED PHOTOMULTIPLIER TUBE SIGNALS, 

0. 25xl0 -6 m PARTICLES, OPTICAL 
EFFICIENCY 0.3 (PULSES LOCATED AT SAME 

TIMES AS IN FIGURE 17) 48 

21 LOW-PASS FILTERED SIGNALS FROM FIGURE 21 49 


-iii- 



LIST OF TABLES 


Table 


1 TABULATION OF MEAN VELOCITY, DATA RATE 
AND APPARENT TURBULENCE INTENSITY 
CURVE FIT 

2 SIMULATION PARAMETERS 



1.0 


INTRODUCTION 


A laser transit anemometer (LTA) measures a two-dimensional 
vector velocity by measuring the transit time of scattering particles 
between two focused and parallel laser beams,. Spectron Development 
Laboratories, Inc. (SDL) has been active in the research and develop- 
ment of LTA systems. The objectives of this research contract were: 
(1) the determination of the concentration levels and light scatter- 
ing efficiencies of naturally occurring, submicron particles in the 
NASA/ Ames Unitary Wind Tunnel Plan and (2) the evaluation based on 
these measured data of a laser transit anemometer with digital corre- 
lation processing for nonintrusive- .velocity measurement in this 
facility. The criteria for the evaluation were to have been the 
speed at which point velocity measurements could be realized with 
this technique (as determined from computer simulations) for given 
accuracy requirements. 

Prior to the NASA contract award, SDL constructed a develop- 
mental prototype LTA system under contract to the Arnold Engineering 
Development Center. The SDL LTA .system was subsequently used to make 
data rate measurements in the NASA Ames ^6-foot x 6-^foot transonic 
wind tunnel during the middle of October., 1978. Specifically, back- 
scatter velocity measurements were' obtained without seeding with mean 
velocity repeatability of 0.15 percent in four successive runs of 
5 seconds each while 'on condition’ at Mach 0.9, Reynolds Number 
3xl0^/foot. Additional measurements were made up to Mach 1.6. 
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In this report we briefly describe the SDL LTA apparatus, the 
measurements at NASA Ames, -the data .collected, and its interpretation 
with SDL computer software and computations, theoretical scaling 
equations for data rate, and. the system design parameters for a 
larger-scale system with a range, of 2 to 3 meters. 
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2.0 


HISTORY AND DESCRIPTION OF SDL LTA SYSTEM 


2.1 Historical Introduction. 

Initial laser velocimeter (LV) experiments in this country 
were made with heterodyne (reference- beam) optical systems using 
frequency tracking electronics.. The dual-beam (fringe) optical sys- 
tems and classical-signal burst-counter electronics have proven to 
be more applicable to high-speed air flows and have been used with 
success. 

During the last several years-, LV systems have been exten- 
sively developed and applied .to practical supersonic and transonic 

flow measurements. These measurements' have been made in the invis- 

1 2 

cid region and within turbulent boundary layers with separation ’ . 

Particle lag problems are sometimes encountered in aerodynamic 
3 4 

LV measurements * . Unseeded (or even filtered) air contains a very 
large number of submicron particles which are sufficiently small to 
follow high accelerations in the flow*, but obtaining adequate back- 
scatter single-particle signal- to-noise ratios for classical burst- 
counter measurements is difficult even with the higher powered argon 
lasers. 

Schodl and others elsewhere have been developing the rather 
old idea of a two-spot system in which the probe volume consists of a 

* W. T. Mayo has reported measurements of ' large numbers of particles 
in the 0.2 - 0.8 micrometer diameter range in the laboratory air 
5 

at the AEDC 5 . 
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pair of parallel highly focused laser beams and a matched receiver 

instead of the now common dual -beam fringes. A. E. Smart reported 

measurements with a different type of two-spot system at Cambridge in 
12 

1977 . The advantages include much less total probe volume and much 

higher focused laser beam intensity; and both of these qualities 
enhances sensitivity to small, particles while decreasing the probabil- 
ity of intercepting the more sparsely distributed large particles. In 
addition, the higher .resolution spatial filtering provides much better 
background light- (flare) rejection when measurements near bodies are 
desired. 

Photon correlation has been known to be ideally more sensitive 
than classical detection for several years. The techniques and theory 
for fringe LV systems are reviewed in. References 6-10. Unfortunately , 
the available correlator electronic speed has not really been adequate* 
to be compatible X'tLth transonic speeds and the simultaneous demands 
for small probe volumes. An LTA system has a great advantage here. 

It achieves lower electronic bandwidth requirements while increasing 
the small-particle SNR and isolation. The reduced bandwidth makes 
the use of a 50 nsec Malvern, correlator feasible for supersonic 
velocities. Thus, the combination of two-spot techniques and digital 
correlation offers exceptional improvements in small-particle sensi- 
tivity and is practical at supersonic velocities. 


* John Abbiss has made supersonic measurements at the R.A.E. in 
England but with large fringe spacing appropriate for the 50 nsec 
resolution of the Malvern Instruments correlator^. 
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A. E. Smart, formerly of Rolls-Royce, Derby, England, has 

taken advantage of these compatible techniques (two-spot and digital 

11-12 

correlation) and extended them in some rather clever ways . He 

designed, constructed, and used a system employing a small helium 
neon laser to make backscatter measurements in both subsonic flows 
and supersonic flows. The details of the Rolls-Royce instrument, 
which is capable of mapping o.ut the .two-dimensional u-v velocity 
probability function for statistically stationary flows, are propri- 
etary to Rolls-Royce and have not been released. However, Dr. Smart 
has joined SDL and assisted W. T. Mayo in the development of a new 
argon-laser transit anemometer system .under Contract F40600-78-C-0002 
for the U.S. Air Force Arnold Engineering Development Centers 

In previous NASA, AEDC,, and ARPA sponsored research, W. T. 

Mayo, Jr. has developed extensive laser velocimeter computer simula- 

15-17 

tion and design software . , photon correlation interpretation 

18 19-21 

software , and specialized photon correlation hardware 

During the period March 1978. through August 1978, SDL has 

developed a completely new laser transit anemometer, system. The SDL 

LTA system utilizes a unique f/,4 mirror-dove rotator/.derotator prism 

for aberration-free spot pair, rotation. The system includes newly 

developed pulse discriminators which operate to estimate the pulse 

center independently of the signal amplitude from classical signals 

down to .photon resolved signals. New fiber opti-c techniques are 

utilized to improve receiver efficiency and probe volume depth 


restriction. The system rotation optics, the discriminators, and 
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the correlator are all under microcomputer control; and the micro- 
computer displays and provides first-order data interpretation to the 
user in semi-real time while storing all raw data for later detailed 
analysis by a larger computer. Figures 1, 2, and 3 illustrate the 
new SDL LTA optical head. The details of the system design are being 
reported separately in the AEDC Contract F40600-78-C-0002 final report 
(in preparation) , but a system overview is presented in the foiloxring 
section which xrill be useful in understanding what follows. 

2.2 Description of System Used for NASA Ames Experiment 

The SDL LTA system,, .partly shown in Figures 1-3, was transported 
to NASA Ames on October 11, 1978, and set up as shoxra. schematically in 
Figure 4. A 50 ns Malvern correlator (48 store) xjas used for. the high- 
speed delay processing of the SDL pulse discriminator, output signals. 
The correlator was used in both the 'photon 1 and the 'pulse' correla- 
tion modes and all the indistinguishable states in between as discussed 

further belox^. The laser was a hexel Model 75-2 operating at about 
-3 -3 

135x10 W at 514.5 nm with 100x10 W in the transmitted beams. 

The parameters of the SDL optical head were as follows. The 

spot separation was measured xrith a calibrated microscope objective 

as 376x10 m with 1/e spot diameters estimated at 1/25 of the spot 

_6 

separation, i.e., 15x10 m. The 'throx-r' or range from the front face 

-3 

of the optical head to the probe volume center x?as 473x10 m. The 

_ q 

transmitter/ receiver lens pair were 40x10 m in radius. The outer 
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Figure 1. Overview SDL Laser Transit Anemometer (LTA) Optical Head. 
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LTA Optical Components 



Figure 2. Schematic Diagram of Optical Components in SDL Optical Head. 



Figure 3. SDL LTA Electronic Components. 




Figure 4. Block Diagram of SDL LTA System for Ames Transonic Measurements 
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annular portion of the lens from radius 20x10 to 40x10 m radius 

-3 

was used for the receiver, while the inner area from 0 to 20x10 m 
radius was used for the transitted beams. The transmitter/receiver 
rotator prism assembly has digital- to-analog DC drive with DC feed- 
back sensing, and a separate linear potentiometer and analog-to-. 
digital readout. The precision of this angular control system is 
0.1 degree. 

The co-annular backscattered radiation is imaged around a 
small transmitter beam turning mirror to two pinhole aperture stops 
which transmit the desired signal .radiation and some diffuse back- 
ground radiation to a microscope objective which magnifies the images 
of the transmitted beams illuminating scatterers onto the ends of two 
fiber optics assemblies. The purpose of the fiber optics is to sepa- 
rate the light from the two spots to two separate photomultiplier 
tube detectors, one for each spot*. The use of two detectors pro- 
vides completely unambiguous direction sensing capability and also 
avoids a zero delay correlation which would interfere with data 
processing in high turbulence conditions. Perhaps the most subtle 
advantage of two detectors is that the effective background light is 
one-half for each detector over .what would occur is both spots were 
imaged on the same tube. Thus, the undesired background for the 
cross correlation of two detector outputs is one— fourth that of the 
autocorrelation of one detector output with both signals. 


* For these experiments EMI 9816 and an EMI D305 tubes were used. 
Both tubes are 2" diameter, 14 stage, S-20 tubes of high quality. 
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The pulse detection discriminators include several overlapping 


ranges. Range 1 is optimized for pulses with duration 25-100 nsec; 
Range 2. for 50-200 nsec; Range 3 for 100-400 nsec, and so forth for 
seven ranges. Each range includes separate optimum detection filters, 
which are flat in response across the range so as not to bias the 
detection probabilities towards . higher or lower velocities, and pulse 
center detection filters which determine the center of the pulse 
independently of signal amplitude. 

for signal pulses which contain only one, two, or a small num- 
ber of single photoelectron .pulses, the pulse center detection filter 
measures a crude weighted average pulse center estimation. This means 
that pulse center jitter errors -can occur up to the transit time dura- 
tion as the number of detected, photons goes down to one. This implies 
a distinct precision advantage of working with three or more photo- 
electron pulses per scatterer. transit both to avoid processing so many 
single photoelectron signals arising .purely from background light and 
to reduce the pulse-center- estimation error. In order to use semi- 
classical and. classical pulse techniques^ one raises the detection 
threshold so that single photoelectron pulses are detected less often. 
However, due to the statistical distribution of single photoelectron 
pulse heights from the PMT electron, multiplier , the transition from 
single to multiple photoelectron events is not distinct, as the data 
discussed later shows. 

The discriminator detection range was both manually and micro- 
computer selectable from the .remote location where the operator 
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station was located. The discriminator level and the PMT high-voltage 
were both manually remote controllable., The PMT's were protected 
from maximum anode current overload by high-voltage relays following 
an anode DC current monitor. These relays were remotely resetable 
(not shown in figure) . The prism rotator also has remote manual,- 
control. Thus', using the correlator on manual with manual rotator 
and discriminator control, setup and checkout operations can be per- 
formed independently of the .microcomputer . When the selector switches 
are set to .microcomputer control,, the control of the correlator start, 
data acquisition, and rotator control, become automatic. 



3.0 


MEASUREMENTS AT NASA AMES 


3.1 Pre-Tunnel Procedures 

The system illustrated in Figures 1-4 was set up by W. T. Mayo 
and A. E. Smart with the optical head located on a table pointing into 
a two— inch thick schlieren quality window upstream of the primary test 
section of the 6' x 6' transonic wind tunnel. It was determined expe- 
rimentally, that less background counts from single-photoelectron sig- 
nals were obtained (~800',000 per channel) when the optical axis was 
tilted slightly to avoid the direct window reflection back into the 
transceiver lens*; it was noted, however, that a perpendicular arrange- 
ment allowed operation with approximately 4x10^ background counts per 
channel. These figures were measured. prior to a change in the optics 
which preceded the measurements of October 17, and are thus only quali- 
tatively useful. All of the final signal data was acquired on Range 2 
(smaller effective gain by 2x) with the PMT voltages reduced to 1900 v 

and 1960 v for the two tubes, with thresholds from 50 mv to 400 mv. 

—1 

The resulting input pulse count rates were less than 200,000/s in 
most cases due to the higher relative detection threshold. 

One of the difficulties of the experiment was that the tunnel 
was inoperative for most of the, first three days after - the equipment 
was set up. A brief run just after setup indicated that data was being 


* These were the final rates on October 17 and are large due to lens 
aberrations. Lower values are expected in the future. 
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obtained at- transonic speeds * but the sensitivity was marginal. This 
fact and the down time provided incentive for a small modification 
of the system optics wherein the microscope objective which images 
the spots onto the fiber optics was reduced in power. Earlier 
laboratory measurements had indicated that the transceiver lenses 
were nearly diffraction limited over the transmitter central position 
but that the effective blur circle of the receiver annulus was ten 
times larger in diameter. Since the design had assumed that the 
return primary image was the matched size, the actual scatter images 
on the ends of the fiber optics were too large, even though some 
reduction had already been accomplished. The additional change of 
imaging objective reduced the image to* only about three times the 
diameter of the fiber. This increased, the sensitivity noticeably. 

To further quantify the. degree .of difficulty arising from the 
loss between the input and output of the fibers, an experiment was 
conducted in which a small NASA-owned low-gain photomultiplier tube 
was used to measure separately the light imaged from a finely diffused 
metal surface in the probe volume center before entering and after 
leaving the fiber optics assemblies. It was difficult to perform this 
measurement precisely due to the necessity of using a dentist mirror 
to get the light that would have imaged into the fibers to image into 
the pinhole of the test PMT.. Also, there was indication that the tube 
was saturating slightly, even at- the lowest available laser power 
setting. In spite of these considerations, the results showed that 
only about l/10th of the available -light was reaching the PMT's. 


- 15 - 



This is consistent with our visual - observation that the scatter image 
appeared. nearly three times ; wider -than the . end' of -the fibers.- We 
have performed laboratory experiments which indicate that the short 
focal length - internal lens has most of the aberrations. This will be 
corrected in future systems. 


Tunnel Meausrements: 


1,, October 16, 1978 


When the tunnel was again operated, it was swing shift on 

October 16, 197,8. The tunnel was operated at a variety of Mach 

6 6 

numbers .from 0.6 to 1.6 at both Reynolds Number 1.5x10 and 3.0x10 /ft. 
The system appeared to perform very well. A variety -of system param- 
eters were adjusted and the .following general selections for each run 
were made: 

Discriminator Range: 2 

Sample Increment: 50 ns 

O 

No. of Samples: 10 

Total Time 5 s 

HV PMT A: 1900 v 

HV PMT B: 1980 v 


Discriminator Thresholds: 50, 70, 100, 140, 200, 280, 400' mv. 

It was observed that easily- discernible signal correlation peaks could 

be obtained with 0.5 s and even 50 ms runs during subsonic runs with 

Reynolds Number = 3x10^/ ft. During the initial run, the Mach number 

£ 

was uprto. Mach 1.6,, Reynolds Number = 1.5x10 /ft, and large amounts 
of condensation occurred. Visibility into the lighted test section 
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downstream was poor. The P.MT protect relay tripped, and the high 
voltage was reduced to 1800 v and then 1600 v. Good signal correla- 
tions were obtained at 1600 v with threshold at 280 mv. 

During the Group 1 experiments , the effects of .threshold changes 
were noted and also the effects of angle rotation. Since the probe 
volume was located approximately 20 cm into the free-stream flow, the 
turbulence level was very low and the angular increments were quite 
small. The data from these runs was presumably recorded, with a few 
auxiliary hand written notes concerning run numbers and threshold 
settings. Unfortunately, examination of the contents of the tape 
recordings the next morning (October 17., 1978) revealed that a.minor 
software error had been made in final changes just prior to transport- 
ing the equipment which not only was not detected in the setup pro- 
cedures but also caused all the data from -the- Group 1 experiments to 
not be saved. Thus, during the day of October 17, 1978, the software 
was modified in preparation, for the experiments during swing shift on 
Tuesday, October 17, 1978. 

3.3 Tunnel Measurements: Group 2, October 17, 1978 

The scheduled tunnel runs for October 17 were to be a sequence 
of many conditions with^only about one minute for each condition. (It 
takes 10-15 minutes to get from one condition to another, but they 
were not going to hold on each condition very long.) However,, it was 
agreed that the condition of Mach 0.9 would be held for 10 minutes 
at each of the Reynolds numbers 1.5 and 3.0xl0^/ft so that SDL 
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could take repeated measurements with the tunnel 'on condition. 1 Data 
sheets were Xeroxed with the genera! parameters previously determined 
already, written in so that the investigators could concentrate on 
writing down the measured count rates on. Channels A and B (using the 
correlator manually) and tape record the correlo grams and other data 
for many different threshold settings at the same condition. Some of 
the data .was taken while the tunnel was changing from one condition 
to another, but several sequences were obtained 'on condition 1 at 
Mach 0.9. at both Reynolds numbers. Generally, the data rates were 
higher for the higher Reynolds number condition. 

3.4 Data Presentation: Group 2 

Upon return to the SDL Costa Mesa facility after the data 
collection on the evening of October 17, 1978, the data tapes were 
transferred to the SDL HP 21MXE System 1000 computer for plotting and 
data .analysis. The raw data. is. .reproduced in the appendix which con- 
tains plots of the correlograms (simple straight line connections 
between points, not curve fits). 

Spectron has software which uses nonlinear regression tech- 
niques to fit data to analytic functions with parameters. Figures 5 
through 11 illustrate the results of first subtracting the mean base- 
line (points outside the peak) and then fitting a 3-parameter Gaussian 
curve (amplitude, mean, and rms width) to the five values near the 
peak. The results indicate dramatically that the determination of the 
mean delay can be much more precise than the 30 ns time delay increment. 
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Figure 7. Curve Fit of Correlogram at M = 0 . 9 , Threshold = 200 mv. 
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Figure 8. Curve Fit of Correlogram at M = 0.9 
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Figure 10. 


Curve Fit of Correlogram at M = 0.9, Threshold = 70 mv. 
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Figure 11. Curve Fit of Correlogram at M - 0.9, Threshold = 50 mv. 


The apparent ms turbulence intensity was typically 2 to 3 percent in 
these runs. These figures are high because they ewe affected by the 
50 ns delay increment of the correlator. The results of applying the 
curve fit software to all of the data runs is summarized in Table 1. 

The fluctuations in apparent turbulence intensity represent the 
effects of tunnel velocity drift over 5-second intervals, photon noise, 
discriminator jitter, etc. The basic width of 4 percent due to the 
50 ns delay increment results- in an apparent rms width of approximately 
2 percent. Figure 12 indicates very small changes of apparent turbu- 
lence intensity over two sets of runs 'on condition' at Mach 0.9, 

£ 

Reynolds number 3x10 /ft. I-t appears that the apparent turbulence 
intensities may be corrected by subtracting the correlator time reso- 
lution effect from the mean square deviation; however, such an analysis 
is beyond the scope of this report. The purchase of a faster correlator 
would be of assistance in making very low turbulence measurements. 

The 'background' presented in Table 1 is the mean value of the 
flat background level in the correlograms (5 seconds of data) . The 
'data rate' is the sum of the correlogram stores near the peak minus 
the background level, the quantity divided by 5 seconds. Thus, the 
sum of the correlogram values above the background in Figures 5-11 
are obtained by multiplying the 'data rate' by 5 seconds for the 
appropriate run number. 

Attention is called to the sequence of runs 75 through 78. 

These four successive runs show a repeatability of the mean flow 
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8 
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47 
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6 
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45 
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33 
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46 
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8 
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6 
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48 
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7 

4B 
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49 

. I488E+88 

3B 
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59 
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6 

0 
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Table 1. Tabulation of Mean Velocity, Data Rate, and Apparent 


Turbulence Intensity Curve Fit. 
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78 

. 7888ESI 

182 

3912 

OVERFLOW 

71 

. 5888ESI 

382 

5422 

OVERFLOW 

72 

. 5888ESI 

218 

5489 

296.342 

.8213 

•7*1 
i" -u* 

. 7898ESI 

■?*Q 

3868 

298.228 

. 8249 

74 

l . 8888ESI 

188 

9 ! 

J uv 

228 . 726 

. 8287 


. I4Q0E+00 

146 

626 

296.848 

.8214 


. 2888E+88 

182 

58 

296.991 

. 8222 


. 2888E+88 

75 

8 

296.861 

. 8285 • 


. 4888ES8 

55 

8 

•pCfT 

t-m . *«•«■ Cm. 

.8224 

79 

. 4888E+88 

119 

8 

/ 29 . 929 

.8168 

89 

. 2828ES8 

176 

! 

OVERFLOW 

81 

. 2988E+88 

189 

57 

198.232 

.8177 
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58 

1 
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92 
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43 

8 
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* OVERFLOW indicates 

that the 

curve fit program 

did not converge. 

In 

come 

cases, the raw 

correlogram in the Appendix shows why in an 

obvious 


maimer. A few cases merit further consideration of the choice of 


initial parameter selection for the software. 


Table 1. Tabulation ;of Mean Velocity, Data Rate, and Apparent 
Turbulence Intensity Curve Fit (Concluded) . 


- 28 - 



RELfiTIVE RMS TURBULENCE INTENSITY 

^ .000 Q.Q1Q 0.020 0.030 0.040 0.050 


.000 


0 X = 65-69 
+ X = 72-78 


+ 


* 


+ 


0 


0 

+ 


+ 

0 


0 

+ 


+ 

0 


1 j 1 1 1 1 i 1 

50 .000 100 .000 ISO. 000 200.000 250.000 300 . 000 350.000 400.000 

THRESHOLD { MV J 


Figure 12 . Apparent Turbulence Intensity on Condition 
at~M~= 0.9, Re.- No. = 3xl0 6 /ft. 


velocity measurement of better than 0.15 percent even though differ- 
ent discriminator thresholds' were used. Other measurements while 
nominally * on condition* show- slightly less stability, but one cannot 
exclude that the tunnel velocity may have been drifting within small 
deviations also . 

Before concluding this . section, we should point out the com- 
parison of our measured mean, velocities with those obtained by the 
tunnel calibration computations. According to the tunnel computer, 
the actual Mach numbers, velocity in -ft/sec, and Reynolds numbers 
during our runs Xi= .65 to 78 were: 


Mach No. V (ft.s Rn/L 


0.899 . 958.27. 2.95 

0.898 957.16 2.94 

0.899 957.83 2.95 


We must take into account that our measurements were made at a dif- 
ferent location than the tunnel probes, and that we independently 
measured our spot separation to probably no better than 1 percent. 
Nevertheless, we have from the tunnel data (approximately the same 
time as runs 75-78 would have been made) : 

• V ‘(ms -1 ) = 957.83 x 0.3048 m/ft = 291.9 m/s 
in contrast with our measured value -of 296.8, i.e,, a difference of 
1.6 percent. 
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4.0 


DATA ANALYSIS AND INTERPRETATION 


Some of the original computational objectives of the proposed 
work have been surpassed by the experimental evidence. The primary 
case of this was our original intention to use single-beam scattering 
amplitude/rate data to predict by simulation some of the accuracy • 
limitations of an LTA with correlation. processing. The truth is that we 
could never have completely simulated all of the error sources and/or 
equipment limitations. The experimental mean-flow repeatability of 
0.15 percent over four successive runs with 5 seconds of data each 
proves that excellent measurement precision is obtainable with a prac- 
tical system. The facts that only an 80 mm diameter lens (whose aberra- 
tions cost an additional power loss of a factor of ten which should be 
recoverable) , and a 135 mW laser was used indicates it-- should be pos.-r 
sible to extend the system to a range of 2 meters (a factor of 4) by 
using 6 to 8 inch diameter high-precision optics. 

Despite the feelings of certainty that the experiment has 
provided, we also need certain analytical and computational results to 
better understand the system performance limitations and to assist in 
scaling the design to larger ranges in the most useful manner. In this 
section we first examine the definition of the words 'data rate' and 
how this relates to different detection methods (photon vs. pulse corre- 
lation vs. pulse timing). This requires an ex ami nation of the nature 
of the pulse signals through simulation and computations- of the system 
sensitivity to submicron particles. We finally consider, with simple 
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equations, the manner in which 'data rate' scales with lens diameter, 
focal range, and laser power. The statistical manner in which the 
detectability of signals reduces in the presence of wall flare is of 
great interest, but we must leave such studies until a later time. 

4.1 Definition of Data Rate 

Figure 13 is a reproduction of two oscilloscope photographs 
taken October 16. For these photographs, the signals from the 
Channel A and Channel B phototube monitors (10 ns pulse rise time) 
were summed electronically for display purposes. The maximum classi- 
cal signals go to just over 2 volts where the preamps saturate. 
Careful observation of the original photographs shows a multitude of 
single photoelectron background pulses at the 50 to 150 mv level. 
Thus, as is obvious in the correlation data of Figures 5 through 11, 
when the pulse detection threshold exceeds the majority of the single 
photo electron pulses , the background level of the correlation func- 
tion goes to zero and only the classical signals which produce the 
correlation spike remain. 

Using the above logic in reverse, one can imagine being con- 
vinced that when the correlation function appears to have little or 
no background, then the count rate of discriminator output pulses 
is equal to the data rate. As an example of this, we see in Figure 6 
that the background is nearly zero. From the raw data for run 
X = 77, we found that the total counts into Channels A and B were 
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(a) 


°p I( p<X)^ n PAG E IS 

°° H QUALrr* 



Mach 0.9, Reynolds No. = 3x10 , 2 Second Exposure. 



(b) Mach 0.7, Reynolds No. = 3x10^, 2 Second Exposure. 


Figure 13. Oscilloscope Time Exposures of the Electronic 
Addition of Channels A and B Photomultiplier 
Tube Output Monitors, October 16, 1978. 
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6942 and 8584, respectively; so after dividing by 5 seconds, we deduce 

—1 —1 

that the input pulse rates were 1388 s and 1717 s , respectively. 

However., we find from Table 1 that the data pair rate was only 75 s . 

Could only 5 percent of the particles be passing through both spots? 

This is not reasonable since .the .tunnel turbulence is so low that nearly 

all of the scatterers which cross one spot should cross the other at 

nearly, the same relative location;, and this is evidenced by the gap in 

the baseline of the photograph in Figure 13(a) and many of the other 

photographs taken. (Unfortunately.,- none are available for run 77.) 

And to make the point a little more dramatically, we can consider run 

—1 

76 (see Figure 7) where the .input pulse rates were 15.24 Ks and 
-1 

12.76' Ks , the correlation background is still far less, than the corre- 
lation peak, and the data pair rate was only 102. .s. ^ . 

The observed data indicate, that a. good correlogram peak to back- 
ground ratio with 100 real" data pairs s ^ in the presence of 15 ,000 s ^ 
Poisson random distributed pulses (both large amplitude sin'gle-photo- 
electron background pulses and single-photon, signals from smaller par- 
ticles which failed to produce a photoelectron for the other spot) . It 
is clear from the above examples that data rate must be defined as the 
rate at which correlated pairs of pulses are detected by the pulse dis- 
criminator. circuits, and that the useful data pair rate may be increased 
within limits by decreasing the. threshold to. include many uncorrelated 
pulses which would cause great inefficiency to any classical pulse pair 
timing circuits with any appreciable dead time. 
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4.2 


Data Rate and Background -Level Equations 


For simplicity assume that all .correlated pulse pairs occurred 

at exactly the same time delay. The rate of occurrence, on either 

channel, due to these 'signal"’ pulses is then A g . The height of the 

'signal' part of the correlogram above .the background level is thus 

A T where the sample duration T is given as 
s 

T = NAt 


where N is the total number -of At clock increments. For turbulent flow 
where the signal delay spreads over more than one "delay increment, then 
the height of the signal part of the correlogram is divided by the num- 
ber of delay increments over .which, -the signal delay is spread. 

Now assume that the rate of additional uncorrelated pulses on 
into Channels A and B respectively are A^ and A^. For locations of the 
cross correlogram other than signal, the correlation function is a flat 
level with expected value given by NP(1,1) where P(l~,l) is the probabil- 
ity of both a 1 in -a present sample of. Channel A and also a 1 in a 
delayed sample of Channel B,. i.e.., for the background locations 

P b d,i) - <y, s + v s + Vb )At2 


For cases where A « A & or A^, as in run 76 for example, P^(l,l) rr A^A^. 
To demonstrate the numbers practically, we consider run 76: 


At = 50x10 9 ' s 
N = 10 8 


X = 15 .2 4x10 3 s 1 
a 

^ = 12. 76x10 3 s -1 

X = 102 
s 

Background ^ (15.24xl0 3 ) (12. 76xl0 3 ) (50xl0 _9 ) 2 (10) 8 = 48.6 

This is to be compared with the run 76 correlator mean value obtained 
as Background = 50 (see Table 1) . 

,In this example the background level is reasonably low and the 
experiment -agrees with theory. We see that data rate cannot be predicted 
by measuring the input pulse rates., A' and X^. it must actually be meas- 
ured with a correlator as was .done in. -the experiments. 

4. 3 Data Rate Scaling with Laser Power, and Optical Efficiency 

Figure 14 is a plot of the true data pair rate versus discrim- 
inator threshold setting for .the Mach 0.9 runs, October 17, 1978. The 
data is replotted in log- log form, in Figure 15. The point from run 72 
(50 mv threshold) appears to have been distorted by the single-photon 
limit wherein continued decrease of. threshold produces no additional 
signal or background pulses. However, due to the statistical range of 
single photon pulse amplitudes, -the limit is not sharply defined. It 
is just as likely that the statistical variability of the correlogram 
with high background rates (see run 72 in -Appendix) is a cause of error 
in the estimation of data pair rates for the lower threshold values. 

The log-log plots would indicate that within the classical sig- 
nal regime, the data rate is proportional to the threshold raised to 



* Actual data measured by the monitor stores of the correlator. 
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Figure 14. Data Pair Rate Versus .Threshold for Mach 0.9 Runs, October 17, 1978. 
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Figure 15. Log-Log Plots of Data Rates from Figure 14. 


the minus 0.833 power; or, equivalently, a factor of 10 increase in 
receiver efficiency or laser power, would result in an increase in 
data rate by a factor of 7., provided the! power law. behavior o.f the 
scattering coefficients versus size remains the same .for the smaller 
more plentiful particles and that the power law for number density 
versus size also extrapolated. For practical design purposes we do 
not really require suctTan extrapolation. This. will become obvious 
in the next sections as we look at. particle size sensitivity Simula-' 
tions and ; geometric scalings effects. 

4.4 System Particle Size Sensitivity Calculations 

Figure 16 is a plot of computer .calculations of the Mie scat- 
tering coefficient for spherical particles with index of refraction 

n = 1.4+jO at 180 degrees backscatter angle. The diameter range shown 

—6 *”6 

is for diameters 0.1x10 to 1x10 m. The relative values at 
— 6 ~6 

0.25x10 and 0.5x10 are seen to be representative (not the bottom 
or top of the curve oscillations) and to differ in magnitude by about 
a factor of 10. For many naturally-produced particulate distributions, 
the number density of submicron sized particles behaves as the inverse 
of diameter cubed in the range of interest down to 0.1x10 . If such 

were the case, we would find that an increase of 10 in laser power or 
a decrease of 10 in detection threshold (classical signals assumed, 
not down to single photon limit) would decrease the detectable diameter 
by a. factor' of 2 which would thus increase the available number by a 
factor of 8. We do not consider it coincidental that this agrees 
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reasonably well with the experimental data-rate versus log-log slope 
behavior j.ust described, since our calculations indicate that the par- 
ticles we were observing went down to about 0.5x10 ^ m diameter as we 
now show. 

SDL has extensive laser velocimeter simulation software which 
has been developed over several years' for NASA Langley, the U. S. Navy, 
and the Advanced Research Projects Agency.' This software is described 
briefly in Reference 15-. We •have used, portions of this.' software and set 
the ’ fringe- visibility 1 and beam intersection angle to zero to simulate 
signals from a one-spot system. - With appropriate inputs, we can use 
these results to see what the return signals from each- of the two spots 
would look like. The simulation includes random single photoelectron 
pulse gain statistics for photomultiplier. Por purposes of display, 
we have assumed that four successive, equally spaced scatterers of the 
same size and composition pass through the center of the scattering 
volume. We have also assumed a background single photoelectron rate 
of 5x10 , which is nearly .an. order of magnitude greater than we 
encountered experimentally with a . 0.1-35 watt laser. The fixed system 
parameters for the simulations are shown in Table 2. The variable 
parameters are discussed along, with .the results which follow:. These 
simulations were performed prior to the experimental trip to Ames. 

The simulated range and laser power exceed the experimental parameters 
of 473 mm and 67 mw/spot,in a somewhat compensating manner so that the 
simulations' provide results approximately correct for our experiments. 
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Table 2. Simulation Parameters. 


System 

Wavelength 

Laser Power 

Optical Efficiency 

Detector Quantum Efficiency . 

Transmitted Beam Diameter 

Range 

Collecting Lens Diameter 
• Probe Volume Spot Diameter. 
Single Pulse Duration 


514.5 nm 

0.1 watt/spot 

Varied (0. 045 ,..0. 30) 

0,15 

25 mm 

762 mm 

80 mm 

20.2xl0~ 6 m 
50.5xl0~ 9 s- 


P article 

Velocity 

Index of Refraction 
Diameter 

Background Light 


400 m/s 
1.4-jO 

Varied (0.25', 

0.50x10 ^ m) 

6 

5x10 Photoelectrons/s 
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First let us consider the case where the product of the trans- 
mitter and receiver transmission efficiency is 30 percent and the 
scatterers are 0.5x10 m diameter. The 30 percent ■ figure includes 
a 75 percent transmission figure which has been measured and an esti- 
mated receiver efficiency of just less than 50 .percent. For such .. 
conditions, the signals out of the photomultiplier (10 ns rise time) 
and the signals after low^pass filtering with a Gaussian filter, are 
shown in Figures 17 and 18. Here the signals are so large that the 
background single photoelectron pulses are negligible in magnitude 
and the 'noise' on the classical pulses is negligible in determining 
the pulse center. 

In practice, we observed that about 90 percent of the received 
-light was not focused into the receiver fibers. We had anticipated 
some such loss, and had run a .simulation with an additional receiver 
factor of 0.15 (which multiplied by 0.3 gives 0.045). The results of 
that, simulation are shown in Figures 19 and 20 where we see that we 
are reaching the limit of classical threshold detection, in the presence 
of single photoelectron background pulses with statistically varying 
(Rayleigh distribution assumed) amplitudes. The location of the four 
signal pulses is the same as in Figures 17 and 18; the other pulses are 
the background single photoelectron pulses. 

In Figures 21 and 22 we show the simulation of 0.25x10 m 
diameter particles with the originally assumed 0.3 optical efficiency 
factor. This is .what we would expect with correction of the spherical 
aberrations which limited the light focused into the fibers in our 
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Figure 17. Simulated Photomultiplier Signals, 0.5x10 m Particles 
Optical Efficiency 0.3. 
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experiments. Now we see that the limit of classical detection is 
approximately a quarter of a micron diameter with the other system 
parameters chosen. This is the limit of sensitivity we practically 
expect once the objective lenses are .replaced. Since most natural 
air contains an abundance of particles at the quarter micron size, 
we expect no problems without seeding for the standard 0.5 m throw 
LTA system. In the next section we consider geometric scaling effects 
for extension of the throw to longer ranges. 

4.5 Geometric Scaling Laws 

Some of the NASA Ames wind tunnel facilities require velocity 
measurements at a range of 2 to 3 meters. These ranges are excessive 
for a system using 80 mm diameter optics if the excellent flare rejec- 
tion capabilities possible with a transit anemometer system are to be 
retained. On the other hand, large diffraction-limited optics with 
large angular aperture become prohibitively difficult to manage. The 
thickness of the windows can become a major source of aberrations. 

Final selection must depend on a detailed design tradeoff, but restrict- 
ing-- the lens diameters to 15 cm has many advantages in compactness of 
packaging and ease of use of the system. In this section we briefly 
explore the effects of geometric size and scale changes on the signals. 

First, if the system is scaled in size while maintaining the 
same f /numbers, then the performance would be unchanged to first order. 
Thus, by doubling the lens diameters and focal lengths, we should 
obtain the same performance directly at a 1 meter range as we obtained 
experimentally at 0.5 m range. 
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Now we must consider the effect of changing the output lens 

focal length (and hence f /number.),* The considerations are identical 

to- those of a fringe mode system. First, if F is the output focal 

2 

length, then the probe volume length increases as F and the beam waist 

as F,, so that the intercept area illuminated and- visible to the 

3 2 

receiver increases as F . . The peak intensity decreases as l/F due to 

2 

the decrease in illumination intensity,, and as an additional l/F due 
to the collector solid angle; .However, the detectability of a pulse 
depends on the number of photo electrons during the spot- transit time. 
The spot transit time is increased, in proportion to F. Therefore, 

3 

the number .of photoelectrons per pulse only decreases as l/F with 
range. 

It would be interesting if. the- increase in the rate of particles 
passing through the probe volume canceled the decrease in signal ampli- 
tude, but this doesn’t happen. Experimentally, we have seen that the 
data rate is proportional to threshold raised to the minus 0.833 power. 
Lowering the threshold is equivalent for data, rate of classical signals 
to raising the detected number of photons per scatterer (inverse 
relationship) . We may thus see that the data rate is proportional to 
increases in detected number of photons to the 0.833 power. 

Combining the above results indicates that the data rate behaves 
as 
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where is the normalizing (experimentally used)::focal length and 
P is the laser power normalized by that used experimentally P^. The 
formula is only valid for small excursions of F from F^ and P from 
P^ because the power law behavior, of the Me scattering coefficients 
and the particle size distributions will change at considerably smaller 
(<0.5x10 ^ m) or larger (1x10 ^ m) particles. 

For a conservative engineering design which keeps the small 

particle sensitivity instead of getting the data rate back because the 

large probe volume cross section intercepts more larger particles, we 

recommend an, increase in laser power and/or detection efficiency to 

3 

compensate for the 1/F effect of a longer focus lens. The larger 
probe cross section will then just provide a higher data rate which 
would prove useful comfort margin' for other tunnel conditions and 
higher turbulence. Fortunately , a factor, of approximately 8 to 10 
increase in detection efficiency is expected due to the use of higher 
performance, lenses without the losses,. due to aberration. Thus, any 
extra laser power will be simply for design margin and. higher perform- 
ance in more difficult measurements. 
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5.0 


SUMMARY AMD DESIGN RECOMMENDATIONS 


1. A laser transit anemometer (LTA) has been shown to be 
capable of 0..15 percent precision of backscatter mean-rflow measurements 
with unseeded transonic flow -in the NASA Ames 6' x 6* wind tunnel at a 
range of 48 cm, with--an 80 mm diameter, lens and a 0.135 watt argon 
laser operating at 514.5 nm. 

2. Analysis and simulation shows that purely classical mode 
pulse detection is wasteful of available data which may be extracted 
by pulse correlation techniques; however, in practical background and 
flare conditions with coaxial argon laser, systems, single-photon corre- 
lation detection is also not optimum due .to the statistical variance 
added by the background. The optimum detection is thus- a semi-classical 
detection system with the detection threshold set at a multiple-photon 
level. Correlation detection is required for optimum sensitivity, but 

5 to 10 ns precision is desirable, for the measurement of turbulence 
intensities at values of less than 0.5 percent for transonic flows. 

3. Scaling equations indicate that a superior system with a 
range of 2 to 3 meters is possible using high-performance 15 to 20 cm 
diameter lenses and a 400 mw argon laser. Higher power may be detri- 
mental to system reliability, would incur Glass IV safety precautions, 
and is not recommended. The engineering of such a system is a natural 
extension of the present SDL ,LTA system. Computer simulation indicates 
that the sensitivity of the experimental system and, hence, the scaled 



*“6 

2 meter system would be for particles .less than 0.5x10 m In diameter. 
Operation of a system at 3 meters would degrade the system performance , 
but good data should still be obtained. 
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APPENDIX 


Plots of the Correlo grams from Runs October 17, 1978. 
All Plots Normalized to Peak Value. 
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